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Abstract 

We calculate the equation of state of nuclear matter for arbitrary isospin-asymmetry up to three 
loop order in the free energy density in the framework of in-medium chiral perturbation theory. 
In our approach Itt- and 27r-exchange dynamics with the inclusion of the A-isobar excitation as 
an explicit degree of freedom, corresponding to the long- and intermediate-range correlations, 
are treated explicitly. Few contact terms fixed to reproduce selected known properties of nuclear 
matter encode the short-distance physics. Two-body as well as three-body forces are systematically 
included. We find a critical temperature of about 15 MeV for symmetric nuclear matter. We 
investigate the dependence of the liquid-gas first-order phase transition on isospin-asymmetry. In 
the same chiral framework we calculate the chiral condensate of isospin-symmetric nuclear matter 
at finite temperatures. The contribution of the A-isobar excitation is essential for stabilizing the 
condensate. As a result, we find no indication of a chiral phase transition for densities p < 2pQ and 
at temperatures T < 100 MeV. 

1 Introduction 

Nuclear matter is an idealized infinite system of interacting protons and neutrons. It is the low 
density and low temperature phase of QCD. In this range of densities and temperatures QCD theory 
is non-perturbative and the active degrees of freedom are no longer the quarks and the gluons, but 
color singlet states, the hadrons. For the description of nuclear matter a variety of approaches has been 
developed during the years, being an important issue for its applications to heavy-ion collisions pQ and 
astrophysics [2]. A powerful tool for this task is represented by Chiral Perturbation Theory, an effective 
field theory that incorporates spontaneous and explicit chiral symmetry breaking. 

The basic idea of our approach relies on the separation between the long-range correlations and 
the short- distance physics |3l HI [5] . While the long- and intermediate-range dynamics in the nuclear 
medium, described by Itt- and 27r-exchange with the inclusion of the A-isobar excitation as an explicit 
degree of freedom, is treated explicitly, the unresolved short-distance physics is encoded in contact 
terms fixed to reproduce selected known bulk properties of nuclear matter. For example, the contact 
term associated to the equation of state of isospin-symmetric nuclear matter is fixed imposing that 
the energy minimum at T = is —16 MeV. The correct saturation density, po — 0.157fm~^, follows 
as prediction. For generalization to isospin-asymmetric nuclear matter one has to introduce two more 
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(a) Isospin-symmetric nuclear matter. 
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(b) Neutron-rich matter with Xp = 0.1. 



Figure 1: Pressure isotherms as a function of nucleon density for symmetric nuclear matter (a) and 
neutron-rich matter (b). They display a first-order liquid-gas phase transition. The dotted lines at low 
temperature in (a) show the non-physical behaviour of the equation of state in the phase transition 
region. The physical pressure is calculated using the Maxwell construction. The dashed lines delimit 
the boundary of the coexistence region. The dots indicate the critical temperature. 



contact terms that are fixed imposing that the asymmetry energy at the saturation point is A(po) — 34 
MeV [5]. Many other ground state and single-particle properties and nuclear thermodynamics emerge 
as predictions. 

Explicit Ivr- and 27r-exchange dynamics are so far treated up to three-loop order in the free energy 
density. Consider symmetric nuclear matter first. The basic ingredient to perform calculations is the 
in-medium nucleon propagator. At T = 0, the nucleon propagator in momentum space has the following 
representation [3]: 



Sn{p) = {^+Mn) 



p2 - Ml + ie 



(1) 



where kp is the Fermi momentum of nucleons and = 939 MeV is the nucleon mass. The first term is 
the known free nucleon propagator, the second term is called medium insertion and takes into account 
the presence of a filled Fermi sea of nucleons by means of the theta step function. The calculation 
is then organized according to the number of medium insertions in the diagrams. The relevant many 
body dynamics come from diagrams with at least two medium insertions. Two-body forces arise from 
diagrams with two medium insertions; diagrams with three medium insertions generate the important 
three body forces. Convergence in the series of medium insertions is realized for ^^7^ -C ~ 1 GeV as 
long as four- nucleon forces are not relevant. 

The free energy is written as a sum of convolution integrals P] : 
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Figure 2: Energy per particle as a function of nucleon density for different proton fractions Xp = Z/A 
at T = 0. The dashed hne shows the trajectory of the saturation point as Xp varies. For Xp < 0.12 the 
energy is always positive. 



where F{p,T) is the free energy per nucleon, /Ci, IC2, /C3 are respectively one-body, two-body and 
three-body kernels, and 



d{p) P 
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is the "one-body" chemical potential defined through the relation 

00 

p = 4 j dppd{p) . (4) 


A{p,T) is the so-called anomalous contribution and vanishes at T = [3]. 

The extension of this chiral three loop calculation to isospin-asymmetric nuclear matter requires 
to take into account the different Fermi seas of protons and neutrons. Consequently we distinguish 
in Eq. ([T]) between proton medium insertion and neutron medium insertion, as well as in Eq. ^ 
between proton and neutron distribution functions. Each diagram now involves the sum of all possible 
combinations of proton and neutron medium insertions with their specific isospin factors. 

In section |2] we present and discuss the equation of state of nuclear matter for different temperatures, 
featuring the first-order liquid-gas phase transition, and the corresponding phase diagram. In section [3] 
we use the scheme above exposed to calculate the chiral condensate of symmetric nuclear matter. 



2 Phase diagram of nuclear matter 

In Fig. ([T| we show the equation of state of nuclear matter as a function of nucleon density p for a 
series of temperatures running from up to 25 MeV. The plot (a) displays the pressure isotherms for 
isospin-symmetric nuclear matter, while the plot (b) is indicative for neutron-rich matter and illustrates 
the pressure isotherms of nuclear matter with proton fraction Xp = 0.1 {xp = Z/A). The qualitative 
behaviour of these curves is reminiscent of the van der Waals equation of state with its characteristic 
liquid-gas first-order phase transition. The dotted lines in Fig. (a) indicate the non-physical behaviour 
of the equation of state in the liquid-gas coexistence region. This part of the curve is replaced by 
the physical one (solid line) using the Maxwell construction. The dashed line sets the boundary of the 
liquid-gas coexistence region, which terminates at the critical point indicated by the dot. For symmetric 
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nuclear matter we find as critical temperature ~ 15.1 MeV. With increasing isospin-asymmetry, 
energy and pressure of nuclear matter at a given density grow in comparison to the symmetric case. At 
the same time the coexistence region shrinks. In plot (b) one can observe that the critical temperature 
has diminished to about 6.3 MeV. 

In Fig. (|2| we study the dependence of the saturation point on the isospin-asymmetry. The saturation 
point is defined as the energy minimum at T = 0. The plot shows the energy per particle of nuclear 
matter at zero temperature for decreasing proton fraction. The saturation point shifts toward smaller 
nucleon densities. The corresponding binding energy reduces with continuity, as shown by the dashed 
line, and vanishes at Xp ^ 0.12. Below this value neutron-rich matter is unbound for any density and 
temperature. 
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Figure 3: T — p phase diagram of nuclear matter for different proton fractions. The dashed line shows 
the evolution of the critical point. 

The thermodynamic information is well summarized by the T — p (temperature versus nucleon 
density) phase diagram in Fig. (|3|. It visualizes clearly how the liquid-gas coexistence region shrinks 
with increasing isospin-asymmetry and disappears at Xp ~ 0.05. Neutron-rich matter with Xp < 0.05 is 
an interacting Fermi gas and cannot undergo any phase transition. 



3 Chiral condensate 

The chiral condensate (gg) is the order parameter of the chiral phase transition. The non-vanishing 
of the condensate is connected to the spontaneous breaking of the chiral symmetry of the QCD la- 
grangian with light quarks. It diminishes with increasing temperature and density. The study of the 
thermodynamics of the chiral condensate is a key issue in order to locate the boundary of the chiral 
phase transition in the QCD phase diagram that takes place when the condensate vanishes. The melting 
of the condensate at high temperature and density restores the chiral symmetry and determines the 
crossover from the Nambu-Goldstone phase to the Wigner-Weyl realization of the symmetry. 

Our starting point is the Hellman-Feynman theorem, which allows to relate the in-medium chiral 
condensate to the derivative with respect to the quark mass of the free energy per particle of nuclear 
matter. Using the Gell-Mann-Oakes-Renner relation [mj^fT^y = —rrig {0\qq\0), one can rewrite the 
quark mass derivative as a pion mass derivative, finding the following formula [B]: 

{qq)ip,T) _ p dFip,T) 

(0|gg|0) /2 dml ■ ^' 
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The quantities and (0 | gg | 0) are taken in the chiral hmit {rriq, m.,^ — )■ 0). The term on the left is the 
ratio between the in-medium condensate and the vacuum condensate, F is the free energy per particle 
of nuclear matter. In our model the pion mass appears as an explicit parameter of the free energy 
and it is possible to calculate its pion mass derivative in a direct way. In the one-body kernel /Ci the 
pion mass dependence is implicit via its dependence on the nucleon mass. The operation of derivation 
leads to the appearance of the pion nucleon sigma term = rn^dM / dm^. For our calculation we 
choose the empirical value cjtv = (45 ± 8) MeV [7j. At T = the chiral condensate assumes a simpler 
representation [8]: 

(0|gg|0) fi\ml\ im^ ) dml ]' ^' 

The term in the round bracket is the contribution of the non-interacting Fermi gas. It is essentially 
proportional to the nucleon density and is the leading term of the chiral condensate. The second term 
is the contribution of the interaction driven by the pion exchange dynamics to the condensate {Emt is 
the interaction energy per particle). 
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Figure 4: Ratio of the chiral condensate of symmetric nuclear matter to its vacuum value as a function 
of the nucleon density p for different temperatures up to 100 MeV. The contribution of the thermal 
pions is included. In the liquid-gas coexistence region at low temperatures (T < 15 MeV) the quark 
condensate decreases linearly. 

In Fig. |4] we show the behaviour of the condensate ratio ([s]) as a function of nucleon density for 
different temperatures from up to 100 MeV. The contribution of the thermal pions is also included |9], 
consisting in a further reduction of the temperature dependent condensate. The effect of the thermal 
pions is appreciable at high temperatures, T > 50 MeV. At low temperature ( T < 15 MeV) one can 
recognize the trace of the liquid-gas first-order phase transition on the condensate curves. Using the 
Maxwell construction one finds that the chiral condensate at a given temperature decreases linearly in 
the liquid-gas coexistence region (see curve at T = 0). The net effect of the interactions is to counteract 
the reduction due to the linear density term. A crucial role is played by the A-isobar excitation. Its 
contribution is large in comparison to that of the other terms and stabilizes the condensate. The bending 
of the curves at large densities in Fig. |4] is a consequence of the inclusion of the A-isobar excitation as 
an explicit degree of freedom. Had we taken into account only the pion exchange dynamics without the 
A-isobar excitation, the system would have reached chiral restoration not far from the saturation point. 
With increasing temperature the effect of the interaction weakens and the linear behaviour is practically 
restored. This recovery is the result of a balance between attractive and repulsive correlations and their 
mass dependence. 
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From systematics in the variation of the condensate with nucleon density and temperature presented 
in Fig. |4] no rapid tendency toward chiral symmetry restoration is seen, at least for densities up to 
p < 2 po and temperatures T < 100 MeV. Being far from the chiral phase transition, our chiral 
perturbative approach is reliable in this range of densities and temperatures. 



4 Summary 

In the present review we have presented recent calculations of the equation of state of nuclear matter 
[To] and of the chiral condensate [6] of isospin-symmetric nuclear matter in the framework of in-medium 
chiral perturbation theory up to three-loop order in the free energy density. 

In our approach long- and intermediate-distance correlations driven by pion exchange dynamics 
are treated explicitly, while the unresolved short-distance interaction is taken into account by fixing 
a few contact terms in order to reproduce some selected bulk properties of nuclear matter. Many 
other ground state, single-particle and thermodynamic properties follow as a prediction of the model. 
Correlations in the nuclear medium include two- and three-body forces. The resulting equation of state 
features the characteristic first-order liquid-gas phase transition of nuclear matter, whose dependence 
on isospin-asymmetry is systematically investigated. We find for isospin-symmetric nuclear matter a 
critical temperature of about 15 MeV. 

The chiral condensate of isospin-symmetric nuclear matter is calculated by differentiating the free 
energy per particle with respect to the pion mass, which appears as an explicit parameter of the model. 
Correlations due to A-isobar excitation play the important role of stabilizing the condensate. As a 
result, we find no indication of a chiral phase restoration for densities up to about twice the density of 
normal nuclear matter and at temperatures below 100 MeV. This is the first calculation of the chiral 
condensate at finite density and temperature that systematically includes in-medium two-pion exchange 
interactions. 

The purpose of the present work is to set realistic nuclear physics constraints for ongoing discussions 
of the QCD phase diagram. 
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